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Abstract 

Transport properties of the multicomponent quantum many-body systems 
obeying Haldane's fractional exclusion statistics are studied in one dimension. 
By computing the finite-size spectrum under twisted boundary conditions, we 
explicitly express the conductivity and the conductance in terms of statistical 
interactions. Through this analysis, the effective charge and effective mass for 
collective excitations are determined. We apply the results for 1/r 2 quantum 
systems as well as for correlated electron systems. 
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I. INTRODUCTION 



In low-dimensional quantum systems, excitations are described by quasiparticles carry- 
ing fractional quantum numbers. One of the well known examples is the fractional quantum 
Hall effect (FQHE) ffl, where quasiparticles are classified by the fractional charge and statis- 
tics. In these theories, fractional quantum numbers arise from exchange properties of the 
wavefunction. Recently, Haldane ]2[ proposed a new concept of fractional statistics based 
on the state-counting of many-body systems, which is a generalization of Pauli-exclusion 
principle. We will refer to this as fractional exclusion statistics. Thermodynamic proper- 
ties were already investigated in detail For example, Wu and Bernard Q formulated 
thermodynamic equations, and showed that the statistical interaction is related with the 
two-body phase shift for Bethe-ansatz solvable models. Their method was generalized to 
multicomponent systems |J, and low-energy critical properties were investigated. 

In this paper, we study transport properties of the multicomponent quantum systems 
with exclusion statistics in one dimension. Transport coefficients are closely related with 
wavef unctions, and may usually be calculated in the Green function formalism. Since we 
have only statistics among particles without explicit wavefunctions, it is not trivial to study 
such properties directly from the definition of statistics. We therefore use a trick to avoid 
this difficulty. Namely, by applying the idea of twisted boundary conditions for the finite- 
size spectrum, we calculate the conductivity and the conductance in terms of statistical 
interactions. We then determine the effective charge and the effective mass, and show how 
these quantities are related with exclusion statistics. 

After a brief introduction of exclusion statistics in the next section, we compute the 
finite-size corrections due to the vector potential for multicomponent quantum systems in 
section 3. In section 4, we then obtain the conductivity and the conductance, and determine 
the effective charge and the effective mass. In section 5, we apply the results for several 
interesting quantum systems. Section 6 is devoted to summary and discussions. 
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II. EXCLUSION STATISTICS 



Let us start with a brief introduction for fundamental properties of exclusion statistics 
It is based on counting the change of the dimension of the one-particle Hilbert space 
when a particle is added to the system, which is explicitly formulated as, 
dD a (k a ) , , , 

dN (kl) = ~ 9a ^ k(X ~~ = ~{9ap( k « ~ k 'p) + S a/3 Sk a k' }, (2.1) 

where D a {k a ) and N a (k a ) are, roughly speaking, the numbers of unoccupied (hole) and 
occupied (particle) states specified by the internal quantum numbers a = (1,2, ...,M) and 
corresponding momentum k a . The matrix g a p, which is called statistical interaction, de- 
scribes correlation effects among particles. For more rigorous definition, see Ref. 0. Simple 
cases g a p{k a — k'p) = go~ a po~k a k'. with g = 1 and g = correspond to free fermions and free 
bosons, respectively, and for general fractional value g, we call it ideal fractional exclusion 
statistics. 

Statistical interactions should be independent of N a [[J, and hence eq. (|2~T|) results in 

D a (k a ) = - E 9ap(k a - k' p )N p {k' p ) + D° a (k a ). (2.2) 

We assume that integral constants are given by D° a (k a ) = D°5 a i or D°, which are referred 
to as hierarchical and symmetric bases, respectively. Such bases are originally used for a 
classification of the FQHE |7],[| . Also, in one-dimensional quantum systems, the hierarchical 
basis serves as a natural basis for the Bethe-ansatz solution |5|]. 

In the thermodynamic limit, we introduce the distribution functions for particles and 
holes, 

Pa (k a ) = ^^ } P f{K) = ^§±, (2.3) 

where D° is proportional to the system size L such as D° = L/2tt under periodic boundary 
conditions HJ. The bare charge for each elementary excitation is defined by |J 

D° 

t a = ^. (2.4) 



Consequently, eg. (pT2|) can be written as 

i-OO 

Pa(K) + pi h \k a ) = dk' p g' a0 (k a - k' p )p p {k' p ), (2.5) 

The energy of the system is assumed to take the form |4| : 

£ = rjo = £ dk a e° a (k«)p*(k a ) (2.6) 

with the bare energy function e°(fc). Note that many-body effects due to exclusion statistics 
are incorporated in the distribution function p(k). 

Thermodynamic equations are generally obtained in a set of coupled nonlinear equations 
at finite temperatures.^] Here we restrict ourselves to the zero-temperature case which is 
sufficient for the following calculation. At zero temperature without external fields, M 
species of elementary excitations are specified by the dressed energy function e a {k). The 
"Fermi level" Q a for each excitation is determined by the conditions, e a (k a ) < for \ k a \ < Q a 
and e a (k a ) > for \k a \ > Q a , provided that the energy dispersion may be symmetric around 
the origin k a = 0. Then, eq.Q2.5p reduces to the following integral equations supplemented 
by those for the dressed energy e a , 

dk '/39ap( k » ~ k 'p)Pp( k '/3)i ( 2 - 7 ) 

p ~ Q P 

e a (k a ) = e° a (k a ) - pt a -J2 [ dk 'f39pa( k 'p ~ k <x) e p( k p)- ( 2 - 8 ) 

j -Qp 

The total energy is now expressed by the dressed energy as 



1 Thermodynamics of systems with exclusion statistics can be formulated by the method proposed 
in Ref. 0, and a multicomponent generalization can be found in Ref. @. A key idea is to introduce 
the entropy S = InW with 

(D a + N a -l)l 
11 N a \{D a -l)\ ' 

which plausibly interpolates the boson and fermion cases Note that if we choose g a p for free 

fermions and bosons, p in eq,([2.3|) reduces to the Fermi and Bose distribution functions. 
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e = 2_^e a + ^n c , e a = t a dk a e a (k a ), (2.9) 

a J —Qa 

where \i is the chemical potential and n c is the density of charged particles. These equations 
can describe static properties at zero temperature. 

III. FINITE-SIZE CORRECTIONS DUE TO STATIC VECTOR POTENTIAL 

We now turn to the computation of the conductivity in pure systems without randomness, 
which can be calculated with the response to the vector potential. Consider the ring system 
at T = threaded by the magnetic flux, which gives rise to static vector potential along the 
ring The effect of the vector potential is incorporated in twisted boundary conditions 



n|,[T2|. Therefore, the energy increment quadratically proportional to A can be calculated 
through the analysis of the finite-size spectrum, which directly gives the charge stiffness and 
hence the conductivity. By observing that the basic equations ( |2.7|) and ( |2.8| ) have the same 
structure as the Bethe-ansatz equations, we can apply elegant techniques of the dressed 
charge matrix developed for integrable models ||13|| , and generalize the calculation of the 
conductivity [[TT],[nj to multicomponent cases. 

Before proceeding with the finite-size corrections, we here give a formal solution to ( |2.7| ) 
in the absence of external fields, which is necessary for the following discussions. To this 
end, let us first introduce the functions |T3|JT^| 



K a/3 (k a - fej) = g' a p(k a - tip) - J2 dk'^ika - tyK jP (k% - k ), (3.1) 

/Q-y 

dk'^Z^ik'^g'^k'^ - hp), (3.2) 

where Z a p is called the dressed charge matrix fL3f . By noting the relation 

/Qa 

dk' a K aP (k a - k' p ), (3.3) 
-Qa 

we can formally write down the distribution function in terms of the dressed charge matrix, 

Pa(k a ) =J2 t l3 Z (3a(k a )- (3.4) 



We now evaluate finite-size corrections to the total energy due to the vector potential, 



which is the same as those for twisted boundary conditions [|lT| , |T2|| . First, by observing that 
the effects of the static vector potential A are to shift the momentum by the amount 5 a 
proportionally to A (in unit e), the basic equations should read, 



p a {ka) 
e a(^a) : 



E 



dk' p g' aP (k a - k'p)pp{k' p ) 



r d k'/39'pa(k'p-k a )ep(k'p), 

-Qa+Oft 



Qa+Sa 



dk a e°(k a )p a (k a ) = J2e a + pn c 



(3.5) 
(3.6) 
(3.7) 



where 



Y + 8a 



J a +8 a 



dk^ai^k^j . 



(3.8) 



It is to be noted that the vector potential not only shifts the momentum uniformly, but also 
can rearrange in general the distribution of the momentum via interactions among particles. 
For clarity, we put tilde for rearranged quantities. 

Let us compute the corrections to the total energy ( p.7|) . By differentiating eq.(|3.6|) with 
respect to 5, we have 

de a (k a 



05, 



0. 



5=0 



by using the conditional equation for the "Fermi level", e a (Q c 
again, we have 

d 2 e a (k a 



(3.9) 

0. By differentiating 



-£p(Qfl)Wf3 a (Ql3 ~ k a ) + 9p a (-Qp ~ k a )}Sfo 



8=0 



/Q~t 
dfyg'^ka - k 



i \ 9%(k') 



8=0 



which can be formally solved by the iteration scheme, resulting in 

d 2 e a {k a 



d5pd5 1 



-e'{Qp){K Pa {Qp-k a )+K l 



f3a\ 



k a )}5f3. 



7' 



(3.10) 



8=0 



where prime in e' a stands for the derivative with respect to k a . In the following, we assume 
the relation g a p{k a — k'p) = gp a {k'p — k a ). Using these results, we have 
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8e a 
85, 



8 2 e r 



85,85, 



5=0 



5=0 



25, 



aP 



Qo 



dk a {Kp a (Qp - k a ) + K i 



Pa{ 



k a )} 



Tp^ 



2t Q e'g(Qp)Z a p5 



Pv, 



(3.11) 

(3.12) 
(3.13) 



where Z a p = Z a p(Qp) is the dressed charge matrix defined in Q3.2|) . Therefore, we obtain 

;h 
8 2 e Q 



the finite size correction Ae = e — e due to the small shift of 5 a as, 



1 



85r85~, 



(5=0 



, 7 , 



(3.14) 



aP 

where V^g = v a 5 a p, and u a is the Fermi velocity defined by t> Q = £' a {Q a ) / p a {Qa)- To derive 
the second line, we have used eqs.(|3.4|) and ( |3.13| ). 

The remaining task is to determine the relation between 5 a and A. The result is quite 
simple (see eq. fl3.23Q ), but still nontrivial, which depends on types of external fields. Here 
we briefly depict how to obtain this relation in a bit general way 0]. For this purpose, let 
us first introduce 



i a (k a ) = I dkp a (k), i Q {k a ) = J dkp a (k), 
G a p(k a - k'p) = J dkg' ap (k - k'p). 



When A = 0, eq.i 



is integrated as 

p "-Qp 



/Qf3 
dk,G a p{k a kp)pp(kp). 
-On 



(3.15) 
(3.16) 

(3.17) 



Notice that the first term in the right hand side is regarded as a bare momentum of the 
system, which can be written as p° a = t a k a . In the presence of the vector potential, the 
momentum is shifted as jP a — * jP a — At a . Then eq.( |3.17| ) is modified to 



dkpG a p(k a — kp)pp(kp), 



(3.18) 



where we denote d a = At a . Note that the following calculation can also be applied for other 
types of external fields if we suitably take d a different from At a . Now introduce k a = k a (k a ) 
such that i a {k a ) = i a (k a ), then we find k a (Q a ) = Q a + 5 a and 



p a {k a )dk a = p a (k a )dk a . (3.19) 



By the use of these relations, we can change the integral variables into k a in eq.( |3.18| ). By 
subtracting both sides of eq. fl3.17p from ( J3. 18[ ) and expanding G a p(k a — k'p) up to the first 



order in {k — k), we get 

/Qp 
dk' p g' aP {k a -k' p )Fp{k' p ), (3.20) 
-Q.a 





where F a is defined by |14 



F a (k a ) = (k a - k a )p a (k a ). (3.21) 
The formal solution to eq. ( |3.20| ) is given by 



F a {k a ) = Y J dpZ Pa {k a ). (3.22) 

P 

As stressed above, this formula is valid for arbitrary types of d a . By explicitly substituting 
d a = At a , and comparing ( |3.21| ) and ( |3.22|) by the use of ( |3.4|) , we get the simple relation, 

k a k a = A = (3.23) 

Namely, momentum shifts occur not only for the charge sector but for all indices a. 

Consequently, by combining ( |3.14[ ) and ( |3.23j ), we end up with the final formula for the 
energy increment as, 

Ae = J2t a (ZVZ T ) af3 tfiA 2 . (3.24) 

a/3 

This completes the calculation of finite-size corrections due to the static vector potential. 

IV. TRANSPORT COEFFICIENTS, EFFECTIVE CHARGE AND MASS 

A. conductivity 



According to the formula ( |3.24| ) for the response to the external vector potential, we can 
obtain the charge stiffness as, 
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D c = J2ta(ZVZ T ) al3 tp. (4.1) 
An important point is that D c is directly related with the conductivity through the relation, 

Rea{u) = e 2 D c 5{u) at u = (4.2) 



according to the linear response theory |10,1] 

In many literatures so far, correlation effects on D c have been considered to modify only 



the effective transport mass m* [10-12]. However, the transport mass is not sufficient to 
describe the correlation effects on transport properties, because quasiparticles such as holons 
in one dimension can carry not only the effective mass but also the effective charge. There- 



fore, in place of ordinary interpretation ||iTi[l2| , we propose the following natural expression 



for the conductivity in terms of the effective charge e* and effective mass m*, 

7166*77* 

Rea(w) = -5(u) at u = 0, (4.3) 

m* 



where n c is the density of charged particles. The charge stiffness calculated in eq. (|4.1|) is 
then related with D c such that 

D c ex 4- ( 4 - 4 ) 

m* 

So, the charge stiffness (or conductivity) is not sufficient to derive the effective charge and 
mass, and another quantity is necessary to determine them. We will show that the effective 
charge can be derived from the conductance. 



B. conductance 

It is known that the fractional charge enters in the conductance for the finite system in 
one dimension. Though it is not easy to calculate the conductance without wavef unctions, 
we can compute it by taking into account a universal property of Luttinger liquids, i.e., the 
conductance is solely determined by the correlation exponent for the charge sector. To this 
end, we first define the critical exponent for the charge density correlation function in the 
asymptotic region, 
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(p(z)p(O)) ~ exp(2tk ( ; ] x)x- 2Mf % (4.5) 



where kp = irn c is the "Fermi momentum" for the charge sector, and is usually given by 

(c) 

k F = Mkp in terms of the ordinary Fermi momentum kp- The correlation exponent f c 
is a function of statistical interactions, which is normalized to reproduce f c — 1 for non- 
interacting systems. An important point is that the conductance G c can be determined by 



f c in a universal way [15 



,e 2 



G c = M-U (4.6) 
h 

This formula can be deduced by observing that the conductance is controlled only by the 
charge degrees of freedom.^ Renormalizing e by f c such that G c = Mee* /h, we can naturally 
define the effective charge as 

e* = f c e. (4.7) 

Therefore, the remaining task is to obtain the critical exponent f c in terms of statistical 
interactions. Following a way similar to the last section, we can derive the exponent f c 
through conformal-field-theory analysis of the finite-size spectrum which has been already 
computed in ||. We thus obtain the renormalizarion factor for the fractional charge, 

f c = ^Y.UZZ T Ut^ (4.8) 

in terms of statistical interactions which are implicitly incorporated in the dressed charge 
matrix. By applying the above formulas to eq.( [1.4[ ), we can also extract the effective trans- 
port mass m*, which turns out to be inversely proportional to the velocity v. The expressions 
( [lTD , ( [4.4|) , fl4.7|) and Q4.8|) are the main results in this paper. 



One can see that the effective charge is determined solely by the statistical interactions, 
whereas the effective mass also depends on non-universal quantities such as the velocity. We 
wish to emphasize that the formula for the fractional charge ( |4.7| ) with ( |4.8|) is universal, 
which holds generally for multicomponent Luttinger liquids. 



^Exponent g without randomness in ref. corresponds to f c in this paper. 
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V. APPLICATIONS 



A. Ideal fractional exclusion statistics 

One of the most remarkable applications of exclusion statistics is that for the ideal case, 
in which the statistical interaction is given in a simple form 



g a p(k a - kp) = G a p5(k a - hp). (5.1) 

This model is known to have close relationship to interesting quantum systems such as the 
FQHE, the 1/r 2 systems, etc. 

If we take a bare dispersion as e Q a {k) = e°(k)t a with e(k) = k 2 /2, the ground state 
configuration is Q\ = Q2 = ■■■ = Qm = Q, and the distribution functions are obtained as 

P 

The number of charged particles is 

n c = 2Qv, v = ^t a G~ptp, (5.3) 

a/3 

where the quantity v is related to the compressibility as k c = 4z/ 2 /n c , which is a one- 
dimensional analogue of the filling factor in the FQHE. In the present model, the Fermi 
velocities for each excitation take the same values 

vi = v 2 = ... = v M = v = (5.4) 

As is shown in Appendix, one finds a simple relation between the dressed charge and the 
statistical interaction, 

ZZ T = G- 1 . (5.5) 

Consequently, the renormalization factor ( |4.8| ) for the fractional charge e*/e = f c is expressed 
as 

/c = iE*«( Gr " 1 )^*/>- M) 



M , 

ap 
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It is to be noted that the fractional charge is now explicitly obtained only in terms of the 
statistical parameters G a p for exclusion statistics. Since the charge stiffness is given by 
D c = n c /2, the enhancement factor for the transport mass in this case is derived 



m*/m = f c 



(5.7) 



from which one can see that the enhancement of m* exactly cancels the renormalization of 
charge e*, in accordance with the translational symmetry. 

Let us now discuss more concrete models for ideal exclusion statistics. For an instructive 
example, we consider the statistical-interaction matrix G in hierarchical basis, which is 
derived from the continuum 1/r 2 model [IB| with SU(M) symmetry fIT|]IS |. This model is 
also related to a fundamental series for the hierarchical FQHE []7L B][B0| . The M x M matrix 
for statistical interactions in this case reads, 



G 



2n+l -1 



-1 2-10 



We have from eq. (|5T3| 



0-12-1 
0-12 

M 



(5.8) 



2n + 1 



V 



2Mn + 1 q 



(5.9) 



■. 1 

' 2 

which corresponds to the filling factor in the case of the FQHE. Observing this, we see that 
the matrix G in this model plays a role similar to the flux attachment in Jain's model for 
the FQHE ||; 



G 



XmxT, 



(5.10) 



where xm is the IQH state with the filling factor M, which is attached by 2n flux quanta xi n - 
In fact, one can find the same matrix as Q5.8D in the classification scheme in the corresponding 
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Abelian Chern-Simons theory for the FQHE ggJUg. Therefore, it is seen that the SU(M) 
1/r 2 model has close relationship to the hierarchical FQHE with v = M/(2Mn + 1). 

According to eqs.( |5.6| ) and ( |5.7| ), it turns out that the effective charge and mass are given 
by taking t T = (1, 0, • • ■ , 0), 

(5.11) 



2Mn + l q 
m m 

m* = — = — . 5.12 

2Mn + 1 q v ' 

Note that the expression for the effective charge (|5.11|) is actually in accordance with that 
for the FQHE. In particular, for one component case G = g, we have e*/e = m*/m = 1/g. 
Note that for free systems, i.e., n = 0, we have e*/e = m*/m = 1. 

B. Correlated electron systems 

It is also instructive to apply the results to one- dimensional correlated electron systems. 
In order to fully describe interacting electron systems in the whole energy range, it is neces- 
sary to consider statistical interactions depending on the momentum in a complicated way 



1 22]] . However, if we restrict ourselves to the low-energy conformal limit, we can still use the 
idea of ideal exclusion statistics. In such low-energy region, the critical behavior is described 
by the Luttinger liquid theory, in other words, by c = 1 conformal field theory. In this 



case, we can introduce 2x2 matrix G a p for ideal statistics in (5.1) in terms of Luttinger 



liquid parameters. Since this model has two kinds of elementary excitations, i.e., spinon and 
holon which have two different velocities, v s and v c , we can choose the hierarchical basis 
as a natural one. By analyzing exactly solved models or Tomonaga-Luttinger models, the 
matrix for statistical interaction can be deduced as ||, 



( K P + l _ x \ 



G 



2K P 



(5.13) 



V -1 2 J 

where K p = Zfx/2 is the critical exponent for the 4k f oscillation piece in the density 
correlation function. Note that G n is related with the charge degrees of freedom, G 2 2 to the 
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spin degrees of freedom, and the off-diagonal elements are regarded as mutual statistics. It 
is to be emphasized here that ( |5.13| ) is the universal formula for correlated electron systems. 



According to eq. (|4.7|), the effective charge is given by substituting t T = (1,0), 

e*/e = K p , (5.14) 



which reproduces the known results for the conductance in Luttinger liquids ||15|| . Also, from 
eq. ([4.1|) , the charge stiffness is found to be 

D c = 2v c K p , (5.15) 

from which we obtain the effective mass using eq. Q4.4l ), 

m*/m = Vp/v c . (5.16) 

where Vp is the Fermi velocity for non-interacting electrons. 

For exactly solvable electron models such as the Hubbard model and the supersymmetric 
t-J model, the critical exponent K p and the velocity of holons were calculated exactly as 
functions of the interaction and the electron density |p2B|-f23J. So, we can discuss the effective 



charge and mass for these models. As for the Hubbard model, K p decreases from 1 to 
1/2 as the Coulomb interaction is increased ||25|| , hence the effective charge decreases with 



the increase of the interaction, as stressed by Ogata and Fukuyama |15|| . Near half filling, 
the effective charge always takes e/2 as far as the Coulomb interaction exists. In the case 
of the supersymmetric t-J model |23|,|24|], the effective charge is e/2 near half filling, but 
as the electron density decreases, it continuously increases and reaches the non-interacting 
value e* = e in the low density limit. Also, we can discuss the effective mass for electron 
systems. The results are essentially same as previously discussed [|12]: the effective mass has 
a divergence property near half filling both for the Hubbard model and the supersymmetric 
t- J model, reflecting the metal-insulator transition. 
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VI. SUMMARY AND DISCUSSIONS 



In summary, we have obtained the transport coefficients, the effective charge and mass 
for multicomponent quantum systems obeying fractional exclusion statistics. Their explicit 
relation to the statistical interaction has been derived in eqs. (|4.1| )~ (O) for general systems 



obeying (|2.1|). We have applied the results for the cases with ideal statistics as well as for 
the conformal limit of electron systems. It has been also pointed out that the statistical 
interaction derived from SU(iV) 1/r 2 models are closely related with Jain's construction (or 
the corresponding Chern-Simons theory) for the hierarchical FQHE. 

It is instructive to note that the effective charge in the ideal case is expressed in an 
extremely simple form (5.6) or (|5.11|) in terms of the statistical interaction G a p, implying 



that the fractional charge in the ideal case directly reflects the fractional statistics. In fact, 
we find an alternative way to derive (|5.11| ) using only the definition of fractional exclusion 



statistics. We briefly summarize how to get them intuitively. In the ideal case, the definition 
eq.(|2.1|) reads 



• a 



-G a l at k a = hp. (6.1) 



dDp Q0 

Now imagine the ground state configuration and make a hole to excite the system. The 
above equation implies that if we make a /3-hole, the number of a-particles decreases by 
the amount of G~l. Then, how many charged particles decrease in all? The answer is 
Y^a^aG~l. Now make a hole at a sector a, i.e., t T = (0, 1, 0, ...0) in the hierarchical 
basis. This corresponds to t T = (0, 1, 1, 1) in the symmetric basis, which may create 
M + 1 — a holes of electrons. Therefore, making a hole with unit charge corresponds to 
removing particles with the charge 



6a M+l-a 



J2 G P*> ( 6 - 2 ) 



which in turn defines the effective charge of the excitation. If we adopt G in eq.( |5.8|) , then 
we have 
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e*/e = 1/q, (6.3) 



where q is defined in eq. (|5.9|) . This result coincides with (|5.11|) . 



Finally, another remark is in order. In section III, we have defined the effective charge 
( |4.7|) apart from the trivial degeneracy M. However, it may be possible to include such 
factor in the definition of the charge. In this definition, eqs. (|5.6| ), ( |5 . 1 1| ) and (|5.12|) are 
modified by the factor M, and eq. (|6.2| ) should be replaced by 



K = eEG£ ( 6 - 4 ) 
and therefore, eq.(|0|) is modified as 

e*Je=(M + l-a)/q, (6.5) 

i.e., explicitly, p/q, (p— l)/q, ■ 1/q. This definition for the fractional charge corresponds to 
that used in refs. 
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APPENDIX A: 



In the case of ideal exclusion statistics, there exists a simpler way to derive the formula 
(5.5) without calculating the dressed charge (|3.2| ). We briefly depict this convenient method 



below. Consider the present system without external fields. Then there are two kinds of 
elementary excitation, i.e., the excitation which changes the number of particles, and which 



carries the large momentum. Following techniques of the dressed charge matrix |T^], the 
excitation spectrum was explicitly evaluated in ||, 

Ae = (v/A)n T (ZZ T )- 1 n + vd T (ZZ T )d, 
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where the vectors n and d denote the quantum numbers for excitations, which label the 
change of particle number and the momentum transfer, respectively. Note that the excitation 
which we seek for in fl3.24|) corresponds to the excitation specified by d. From the above 
formula, one can see that these two kinds of excitations are related to each other reflecting 
modular invariance. So, we can easily deduce the excitations for the d-sector once we can 
calculate those for the n-sector. The calculation for the latter excitation is much simpler 
than the former. Let us then calculate the latter by changing Q a — > Q a + AQ a . Then both 
n a and Ae are given by functions of AQ a , and a simple calculation gives As as a function of 
n a such that Ae = (v/4)n T Gn. Comparing these results, we end up with the formula Q5.6|) . 
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